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ABSTRACT 

We present observational evidence for oscillations of magnetic flux density in the 
quiet areas of the Sun. The majority of magnetic fields on the solar surface have 
strengths of the order of or lower than the equipartition field (300-500 G). This results 
in a myriad of magnetic fields whose evolution is largely determined by the turbulent 
plasma motions. When granules evolve they squash the magnetic field lines together or 
pull them apart. Here we report on the periodic deformation of the shapes of features 
in circular polarization observed at high resolution with Sunrise. In particular, we note 
that the area of patches with constant magnetic flux oscillates with time, which implies 
that the apparent magnetic field intensity oscillates in antiphase. The periods associ- 
ated to this oscillatory pattern is compatible with the granular life-time and change 
abruptly, which suggests that these oscillations might not correspond to characteris- 
tic oscillatory modes of magnetic structures, but to the forcing by granular motions. 
In one particular case, we find three patches around the same granule oscillating in 
phase, which means that the spatial coherence of these oscillations can reach 1600 
km. Interestingly, the same kind of oscillatory phenomenon is found also in the upper 
photosphere. 

Subject headings: Sun: oscillations — Sun: atmosphere — Polarization 
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1. The magnetically, dynamically active quiet Sun 



As we enhance the spatial resolution and sensitivity of spectro-polarimetric measure- 
ments, it becomes increasingly evident that magnetic fields permeate the whole quiet Sun 



( Orozco Suarez et al. 



2007 



Martinez Gonzalez et al. 



2008 



Lites et al. 



2008 



Danilovic et al. 



2010|) . The observed Zeeman polarization signals in the qui etest areas of the Sun are highly 



dynamic, having timescales compatible with the granulation ( iLin & Rimmele 



2007 



Danilovic et al 



19991 : 



Harvev et al 



2010|) . These weak magnetic flux concentrations emerge somewhat prefer- 



the so 


lar atmosphere ( 


Martinez Gonzalez et al. 


2007 


Centeno et al. 


2007; 


Orozco Suarez et al. 


2008; 


Martinez Gonzalez & Bellot Rubio 


2009 


; Gomorv et al. 


2010 


; Danilovic et al. 


2010 


). 



In many of these features the field is organized in the form of ^^-shaped loops and survives 
convective motions, the footpoints of these loops being passively a dvected by granular motions 



(Martinez Gonzalez & Bellot Rubio 



20091 : 



Manso Sainz et al. 



2010h . Importantly, some of these 



small-scale loops do reach higher layers, at least 
very quiet photosphere to chromospheric magneti 


the 


ower chromosphere, hence linking the 


sm ( 


Martinez Gonzalez & Bellot Rubio 


2009; 


Martinez Gonzalez et al. 


2010; 


Wiegelmann et al. 


201( 


)). 



In this letter we present another interesting feature of the small-scale solar magnetism. We 
find that both the area and the magnetic flux density of circular polarization patches displays 
quasi-periodic oscillations at constant magnetic flux. 

Oscillations of line-of-sight velocity and brightness of small-scal e magnetic feature s 



in solar net\york, fa c ulae and plage areas ha v e been detected by e.g. 



Kriiger et al. 



(2001) 



Bloomfield et al 



De Moortel et al. 



(l2006h : 



Vecchio et al. 



(12002): 



(120091) : 



De Pontieu 



Jess et al. 



magnetic field oscillations remains ellusive. Recently, 



et al. ( 

zool. 



(2003) 



ites et al 



(1993); 



Centeno et al 



(120091) : 



(20091). However, a clear detection of 



Fujimura & Tsunetal tOO% used HINODE 



spectro-polarimetric data with a spatial resolution of 0.32", showing evidence of magnetic flux 
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oscillations, interpreted by these authors as an indication for magneto-acoustic "sausage" and 
"kink" modes in individual elements of plages and pores. Very little (if anything) is known about 
possible oscillations of magnetic features in solar internetwork areas, with much less magnetic 
flux. Due to the usually low amplitude of the polarization signals, the observations necessary to 
study the temporal behavior of small-scale patches in the internetwork are rather demanding. Here 
we make use of the superb quality of the Sunrise/TMaX data and show that magnetic fields outside 
active regions also oscillate, representing the first evidence for magnetic field oscillations in the 
quiet Sun intemetwork areas. 



2. The IMaX data 



We analyze disk cente r quiet Sun observations obtained with the IMaX i nstrument 



(Martinez Pillet et al. 



2010 



Solanki et al 



20101) onboard the SUNRISE balloon home observatory (jBartholetaL 
201o[) . The IMaX instrument is a Fabry-Perot interferometer with polarimetric 



capabilities at the Fe i line at 5250.2 A. The observational sequence consists of five filtergrams 
taken at +40, +80 and +227 mA from the Fe i 5250.2 A line center. The field of view is 
46.8" X 46.8", with a cutoff frequency of (0.15"-0.18") \ the pixel size being 0.05". We analyze 
two time series of about 23 and 40 min duration with a cadence of about 32 s and a noise level of 
(T = 10"^ and 7 x 10""^ Ic in the circular and linear polarization (Ic being the continuum intensity). 

We analyze the IMaX data before the phase-diversity reconstruction process, due to their 



lower noise level. The Fe i line at 5250.2 A is a forbidden but resonant 



19871 : 



i ne and thus is ext remely 



Lagg et al.ii201Q) . We 



sensitive to the local temperature of the plasma (e.g. Stenflo et al.l 
use the flux density (0 in Mx cm^) derived from the weak field approximation instead of the 
directly observed Stokes V signal. The weak field approximation states that the Stokes V profile is 
proportional to the derivative of the intensity profile and to the magnetic flux density. Taking into 
account the simultaneous information contained in the intensity and the Stokes V profile, we can 



Fig. 1. — Time evolution of a weak circular polarization patch (corresponding to the top right 
panel of Fig. O. The black and white background represents the magnetic flux density computed 
in the weak field approximation, saturated to +20 Mx cm"^. Blue lines represent iso-magnetic flux 
densities of -130,-100,-75,-50,-40,-30,-20,-10, and -5 Mx cm"^. Yellow line is the iso-magnetic 
flux density contour containing a time-constant magnetic flux of -4.5x10^^ Mx. 

decouple the thermodynamical and the mag netic efi'ects. The magnetic flux density is estimated 



following a least-squares minimization (see 
as: 



Martinez Gonzalez et al. 



20111). It can be expressed 



y • — y,- 



which uses all the wavelength samples across the profile (index z) and is therefore more accurate 
than determinations based on single magnetogram measurements. The constant C is defined as 
4.6686 X lO^'^/lpg, Ao being the central wavelength (in A) and g the effective Lande factor of 
the transition. The derivatives of the intensity profile have been computed numerically using the 
observed data. For this kind of IMaX data, this approximation is valid for magnetic fields up 
to ~ 1 kG. The estimated error is cr^ = 4 Mx cm"^. If we assume that the magnetic fields are 
resolved (filling factor is one), Eq. 1 gives the longitudinal magnetic field strength in G. 



3. Oscillations of quiet Sun magnetic fields 



We focus on the temporal evolution of circular polarization patches. Magnetic patches are 
passively buff"eted by granular motions and interact with them, fragmenting, mixing with patches 
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Fig. 2. — Same as Fig. [T]for the top left panel of Fig. |3l Yellow line is the iso-magnetic flux 
density contour containing a time-constant magnetic flux of -5x10^^ Mx. 



of opposite or e qual pq 



resolution (e.g. 



Martin 



arity, and cancelling. This behavior has been also observed at lower spatial 



1988|) . but at higher resolution we observe more polarization patches and 



hence a large probability of interaction. In fact, we estimate that roughly 50% of the circular 
polarization signals suff'ers at least one encounter with a neighboring patch in 5 min. 

The spatial resolution of the IMaX data (~ 0. 15") allows us to observe an interesting dynamic 
property of the weak quiet Sun polarization patches. The signals are found to vary in shape, i.e., 
the area enclosed in a contour containing a constant magnetic flux changes with time. This can be 
seen in Figs. [Hand [21 following the yellow contours. Both the area and the shape clearly change 
with time. 

Although this behavior is common to all magnetic flux density patches (the ones with 
significant spatial coherency), we select some particular examples to study the area fluctuations in 
detail. We study four low flux patches (10^^ - 10^' Mx) that do not interact with another patch 
during the whole sequence or for a long time (at least at our present noise level) and an intense 
patch (10'^ Mx). We follow the time evolution of these patches at constant magnetic flux. In other 
words, the patch is defined as the signal enclosed in the magnetic flux density contour for which 
the magnetic flux is constant with time. At all times, the contour value is above 8 Mx cm"^, which 
is twice the noise level of the magnetogram. The areas are enclosed in contours that are treated as 
polygons, hence, we increase the precision on the computation. Since the magnetic flux density 
has an associated error, the contours (and the areas) have also an uncertainty. We have estimated 
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the errors associated to the areas numerically using a monte -carlo approach. The larger the area 
the smaller the relative error. In our case, the errors range between 1-8%. These uncertainties 
have been propagated as gaussian variables to derived quantities, such as the mean magnetic flux 
density enclosed in an area (see Fig. |5]). 

Figure |3] shows the time-dependence of the area of four selected weak magnetic features that 
show clear periodic fluctuations of the area. The area fluctuates between a factor of 7 and a factor 
of 2. The four selected cases present a large variety of behaviors. From the top left, clockwise, 
the modes are amplified, damped and showing two different periods and different damping or 
amplification factors. The periods range between 4 and 11 minutes, showing that there is not 
a characteristic period for these oscillatory patterns. These periods have been obtained using 
analytical fits to the area oscillation since the duration of these events does not allow a reliable 
fourier analysis. The corresponding anti-phase oscillations of the magnetic flux density have peak 
to peak variations from 12 up to 100 Mx cm~^ around the mean value, which are -24, -67, 48 and 
-34 Mx cm~^ from the top left of Fig. |3l clockwise. 

Figure |5] shows the area oscillations of a more intense patch (with a maximum value of 
270 Mx cm"^) at three different radii from the center of the structure, defined by the amount of 
flux contained within the considered perimeter. As for the weaker flux patches, the oscillation 
contains several periods. In this case, 3 and 5 minutes are the most evident. The area oscillation 
mantains period and phase over different fluxes. The only difference regards its amplitude, which 
increases as we move away from the center of the patch. The fluctuation of the area relative to 
its mean value ranges from 14%-24% from the most central to the most peripheral contour, i. e., 
the areas containing 1x10''' and 1x10'^ Mx, respectively. The magnitude of the fluctuation of the 
magnetic flux density remains constant. This means that the variation of this quantity across the 
structure is very smooth. 

The line of sight velocity, the continuum and the core intensity oscillate too, showing the 
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Fig. 3. — Evolution of the area of four selected magnetic patches. Constant amounts of magnetic 
flux are followed in time, being -4.5x10^^, -1.0x10'"', 5x10'^, and -9x10^'' Mx, from the top left 
panel to the bottom right one. The measured areas are represeted as error bars. The blue lines 
are analytical fits (using a least-squares estimator) to facilitate recognition of the quasi-periodic 
pattem. The time t in these expressions is given in minutes. 

characteristic pattem of the p-modes, with a 5 minutes period. The continuum and core intensity 
fluctuations exhibit a phase shift of 180°. The velocity is also in anti-phase with the continuum 
intensity. After removing the eff"ect of p-modes using a subsonic filter, the speed of sound being 
4 km s"^ Neither the oscillation of the velocity, the continuum nor the core intensity seem to be 
related to the magnetic flux density (and area) fluctuations, as can be seen in Fig. [3l 
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Fig. 4. — Evolution of diverse physical quantities of the four patches of Fig. |3l The black lines 
represent the oscillation of the magnetic flux density with respect to its time average. Blue, green 
and red lines represent the variation of the line of sight velocity, continuum and core intensity with 
respect to its mean value, respectively. Note that these three last quantities have been multiplied 
by a factor for visualization purposes. Dashed lines show the time interval that is aff"ected by 
apodisation. 

Whereas some of the oscillations in area we observe display periods in the range of the 
p-modes, others have distinctly longer periods, even up to 1 1 min, and show an abrupt change of 
the wave period at a given time. Moreover, we find that the area of some magnetic patches located 
close together oscillate in phase. In one particular case, the co-oscillating flux patches surround a 
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granule, the mean distance between them being about 1600 km. 



4. Discussion 



We have shown that the magnetic field in the quiet Sun oscillates with a variety of periods 
compatible with the granular life-time. These oscillations can be strongly damped (top right panel 
of Fig. [3]) or amplified (top left panel of Fig. |3]) with a damping time of ~ 5 - 30 min. We have 
also observed that the period of the oscillation can change abruptly in the course of the time series. 
This would suggest that either the characteristic modes of the structure vary due to changes in 
their physical properties (density, magnetic field orientation, sub-resolution structure,etc.), or that 
they are not characteristic modes of oscillation at all. 

Theoretical stud ies of the characteristic oscillations modes of thin flux tubes (for example 



the classical work by 



Edwin & Roberts 



1983|) suggest a dependence of the oscillation period on 



the structure size a multiplied by the wave number k, in a sense that structures with larger ka 
should oscillate with larger periods. The wave modes excited in flux tubes by convection can, in 
principle, h ave a variety of wave numbe rs k. However, from the theory of the wave excitation (see 



for example iGoldreich & Kumaij(|1990|) ') one may argue that maximum of the excitation happens 
at wavelengths close to the local scale height in the atmosphere, providing that the modes with the 
wavelength about the tube radius ((~ 250 km in our case) will have mo re power. This argument 



(Il983h . 



allows us to fix k and to speculate that, in principle, based on the work of Edwin & Roberts 
we can expect that magnetic structures with larger radius a would oscillate with a larger period. 
Nevertheless, the strength of this dependence is difficult to estimate in our case, as one should 
know precisely the thermodynamic and magnetic parameters of the medium inside and outside 
the magnetic feature. In any case, the time-average size of the magnetic features we study in this 
letter is approximately the same and, thus, the presence of characteristic oscillation modes would 
not explain the scatter in the temporal periods we observe between the different features. It will 
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also be difficult to explain the abrupt change in the periods of the oscillation of a single feature 
and why some oscillations are damped while others are amplified. But more importantly, the 
variations in the area we detect are so strong that, possibly, a non-linear analysis of the variations 
is required, hence, the classical models are not suitable to explain these observations. 

If not being flux tube oscillation modes, the pattern we observe might be the forcing 
due to the evolution of the granules. Since the magnetic fields on the quiet Sun have mainly 



strengths lower than the equipartition field in the photosphere (~ 300 - 500 G; se e e.g 



Khomenko et al.l 



20031 : 



Martfnez Gonzalez et al.,.2008i: iQrozco Suarez et al. 



Lin 



1995 



2007 ). the field lines 



are continuously buffeted by granular flows, being squashed or released everywhere and hence the 
magnetic fields are constantly being amplified or weakened. However, the present observations do 
not allow us to reject either scenario. 

Now, an interesting question arises: do these waves propagate up through the solar 
atmosphere? We predict that, in case these waves propagate across the solar atmosphere, the 
very same phenomena reported in this letter should be observed in the filter-polarimeters onboard 
the Hinode satellite. The different instruments onboard t he Hinode satellite provide valuable 
simultaneous data tracing the solar atmosphere (see e. g. IMartmez Gonzalez & Bellot Rubio 



20091). This data may be interesting to future work devoted to a detailed study of the propagation 
of these waves from the photosphere to the chromosphere. Wave propagation is an efficient means 
of carrying energy between different atmospheric heights and of dissipating it efficiently, mainly, 
through the formation of shocks. Therefore, we should put additional efforts into unveiling the 
role of this newly detected magnetic field oscillation in heating atmospheric layers. 
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suggestions. The Spanish contribution has been funded by the Spanish Ministry of Science 
and Innovation under projects ESP2006-13030-C06, AYA2009-14105-C06 (including European 
FEDER funds), and AYA201 0-1 8029 (Solar Magnetism and Astrophysical Spectropolarimetry). 
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Fig. 5. — Evolution of the area (left) and magnetic flux density (right) of an intense patch (with a 
maximum value of 270 Mx cm"^). The red line contain 1x10'*^ Mx. Note that the oscillation in 
area of the red line has been divided by 4. The black line represent the fluctuation of the area that 
contains a constant magnetic flux of 3x10^' Mx. The blue and green lines are the areas containing 
2x10^^ and 1x10^' Mx, respectively. The larger the flux (area) the further the contour lies from 
the center of the structure, which is defined as the position of the maximum flux density. 
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